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Abstract. We relate the semiclassical asymptotics of the 6j symbols for the 
quantized enveloping algebra U q (s\.2) at q a root of unity (resp. q real positive) 
to the geometry of spherical (resp. hyperbolic) tetrahedra. 



1. Introduction 

Let r > 2 be an integer, q = exjp(jri/r) and Uq(sl 2 ) the quantized enveloping 
algebra for the Lie algebra sl 2 . The category of finite dimensional representa- 
tions of U q (sl 2 ) has a semisimple subquotient J-(U q (sl 2 )) called the fusion cat- 
egory (U Section 3.3]. The isomorphism classes of simple objects in T(U q (s\. 2 )) 
have canonical representatives Vj labelled by half- integers j G [0, (r — 2)/2]. Given 
iii 32 G [0, (r — 2)/2] n Z/2, the tensor product Cg> Vj 2 is isomorphic to the direct 
sum of objects Vj 3 where the sum is over j'3 satisfying the quantum Clebsch-Gordan 
inequalities 

(1) max(ji - j 2 ,j 2 - ii) < 33 < min(ji + j 2) r - 2 - ji - j 2 ) 

and the parity condition ji + j 2 + 33 G Z. Geometrically, the condition Q means 
that there exists a triangle in the unit sphere with edge lengths j a (r — 2)/2n, a = 
1, 2, 3. Generalizations of these inequalities to Lie algebras of higher rank are de- 
scribed in [I], 0,0,123 • 

This paper concerns a generalization of this relationship in a different direc- 
tion, namely from triangles to tetrahedra. The quantum 6j symbol is a function of 
a 6-tuple jab, 1 < a < b < 4, defined as follows: Since the tensor product of simple 
modules is multiplicity-free there is a projective basis for Hom[/ 5 ( S [ 2 )(Vj 14 , Vj 12 <g> 
Vj 23 ®Vj 3i ) associated to each way of parenthesizing, parameterized by half-integers 

{7l2 723 7l3 1 
> are the coefficients in the change 
J34 314, 324 J 

of basis matrix. 6j symbols for q — 1 were introduced as a tool in atomic spec- 
troscopy by Racah , and then studied mathematically by Wigner 021 • 6 j sym- 
bols for U q (si 2 ) were introduced by Kirillov and Reshetikhin |22|. who used them 
to generalize the Jones knot invariant. Turaev and Viro used them to define three- 
manifold invariants 39 , or what physicists call quantum gravity with cosmological 
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constant |JjJ. Roughly speaking these invariants are obtained by combinatorial in- 
tegration of a product of 6j symbols: for the Turaev-Viro invariant, a 6j symbol 
is attached to each tetrahedron in a triangulation, while in the Jones invariant a 
6j symbol is attached to each crossing. Various outstanding conjectures concern 
the asymptotics of these invariants: The Witten conjecture relates the asymptotics 
of the Turaev-Viro invariant (the norm-square of the Reshetikhin-Turaev invari- 
ant) with Chern-Simons invariants of flat 5'0(4) bundles |14) while the volume 
conjecture of Murakami-Murakami relates the asymptotics of the colored Jones 
polynomial with the hyperbolic volume of the knot complement |3l)|. A natural 
approach to these conjectures is to show that non-Euclidean geometry appears in 
the asymptotics of the 6j symbol; we will show that this is indeed the case. 

The connection of the 6j symbols to geometry arises as follows. By a theorem 
of Finkelberg 13 , T{JJ q {si2)) is isomorphic to the tensor category of level r — 2 
representations of the affine Lie algebra s^. The product for the latter category 
uses as its definition the space of genus zero conformal blocks for Wess-Zumino- 
Witten (WZW) conformal field theory. A picture is perhaps the best way of getting 
across the idea of the role the quantum 6j symbols play in WZW: 




The vector space of conformal blocks is isomorphic to the space of holomorphic 
sections of the determinant line bundle on the moduli space of flat SU (2)-bundles, 
see Pauly |31| . The 6j symbol is equal to the Hermitian pairing between two 
holomorphic sections. According to geometric quantization, the sections should be 
interpreted as quantum -mechanical states, and one should try to show that in the 
semiclassical limit the states concentrate to Lagrangian submanifolds. Assuming 
this holds, the leading term in the asymptotics of the pairing is a sum over intersec- 
tion points of the Lagrangians, involving Poisson brackets of the functions defining 
the Lagrangians and the pairing between the sections in the fiber of the determi- 
nant line bundle, see Borthwick-Paul-Uribe [HJ. The particular case of interest is the 
moduli space of flat SU (2) bundles on the four-holed two-sphere, with holonomies 
around the four boundary components fixed by the labels ii2, J23, J34, ji4- The 
Lagrangians consists of flat bundles for which the holonomy around the interme- 
diate circle correspond to J13, J24 respectively, as in PJl. Using the diffcomorphism 
SU(2) — > S 3 , the intersection points correspond to spherical tetrahedra with edge 
lengths l a {, :— Trjab/(r — 2). The Poisson bracket turns out to be the determinant 



In this subject the term quantum is used rather confusingly in two senses: q ^ 1 and H 7^ 0. 
Prom our point of view, — ln(<jt) 2 describes the curvature or, in physics language, the cosmological 
constant. The true quantum parameter is h, which appears in the guise of l/k where quantum 
states are sections of the fe-th tensor power of the determinant bundle. 
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Figure 1. The spherical tetrahedon r 



of the Gram matrix of the tetrahedron, and the phase shift is the area of a holo- 
morphic disk in the moduli space which can be computed using Schlaffi's formula. 
These ideas lead to the conjecture that if r is non-degenerate then 

kju kj 23 kj 13 \ 27rcosO(fc) +tt/4) 

fej'34 



(3) 



kju 

as k — > 00 where r(k) = 

4>{k) 




q—cxp(-iri/r(k)) 

2) + 2, 



r(fc) 3 /a det(cos(U)) 1/4 



^e o6 (fe)/ o6 (fe)-2Vol(r(fe)) 



r(fc) is the spherical tetrahedron with edge lengths 

kjab ~t~ 9 



lab{k) = 27T- 



and 9 a b(k) are its exterior dihedral angles. The formula © generalizes one for 
the case q = 1 given by the physicists Wigner 02| , Ponzano and Regge and 
proved by J. Roberts [32], see also [TT],|2I]. The numerator of the formula was 
conjectured by Mizoguchi and Tada |23|. See also Roberts' discussion 

Our proof of © is not based on the above ideas. Instead, we show (following 
Schulten-Gordon [3^1 and [27j) that both sides satisfy a second order difference 
equation as one label is varied. It follows that each side is a linear combination of 
the two linearly independent solutions to the equation. Taking the Euclidean limit 
and applying Roberts' theorem |34) shows equality of the coefficients. This method, 
while less geometric, has the advantage that it works for non-integral r and gives 
a formula for degenerate tetrahedra with non-degenerate faces. In fact there are 
seven different cases depending on the type of degeneracy of the tetrahedron r; 
these are pictured in Figure 2. The second column, titled Maslov diagram, shows 
the geometry of the relevant Lagrangian submanifolds (or, for the more degenerate 
cases, isotropic submanifolds). The possible degenerations in the table are 

(a) -> (6) - (c) -> ((d) or (e)) - (/) -> (9). 

We refer to case (a) as the non-degenerate case, and case (b) as the tangent case. 
There is an overlap between cases (d) and (e), when the vertices are colinear but 
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Case Tetrahedron Maslov Diagram Power of k 
(a) /X^T /Y/^^\\ - 3 / 2 



(b) 

(c) 

(d) 

(e) 
(f) 




-4/3 

-5/4? 

-1 

-1 
-1/2 



(g) o . 

Figure 2. Seven types of tetrahedra and the corresponding Lagrangians 

I23 = 0; perhaps cases (d) and (e) should not really be considered separate. We 
managed to cover them almost completely. The exception is the case (c) that 
exactly one face is degenerate: although one can compute the asymptotics, we 
did not find a simple formula or manage to construct non-trivial cases where the 
configuration occurs. We also failed to cover the eighth "classically forbidden" 
case that the tetrahedron does not exist, where one expects (and numerically 
experiments show) exponential decay of the 6j symbols. The main theorem (case 
(c) is intentionally missing) follows: 



Theorem 1.0.1. Let r > 2,j 12 ,...,j M E [0, (r - 2)/2] n Z/2 and 
s(k) : { J* a t J23 ) 



kji2 kj 23 kj 13 
kj 3i kj 14 kj 24 J g=exp(Wr(fe)) ' 



(a) If t exists and is non-degenerate, then 

27rcos((/>(fc) +tt/4) 



s(k) 



r(k) 3 / 2 det(cos(U)) 1/4 ' 
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(b) If t exists, is degenerate and all faces are non-degenerate then 

'2\ _1 cos(fc 6 ab j ab ) 



s(k) ~r(fcr^2§3-V/ 3 r M 



3J (AxAzAsAi) 1 /* 

where A a — det(cos(ib c )b !C ^ a ) 1 / 2 and cos(fc ^ dabjab) — ±1 is a s*5« de- 
pending on the sum of the j a t 's for which 6 ab = it. 

(d) If t exists and has exactly one edge length, say l ab vanishing, then 

s(k) ~ (-l) fc ^ + ^ + ^Vr(fc)- 1 (sin(Z ac )sin(/ M ))- 1/2 . 

Note that r need not have any non- degenerate faces. 

(e) // r exists, all faces are degenerate, the vertices lie on a geodesic of length 
at most 7r in order a, 6, c, d and the edge lengths l ac , l b( j ore non-vanishing, 
then 

s(k) ~ (-l) 2fc -'-^r(fc)- 1 (sin(L c )sin(/ M ))- 1 /2. 

Note the edge length l bc may vanish. 

(f ) If t exists, all faces are degenerate, either 

(i) all non-zero edge lengths are equal and less than tx , or 

(ii) one edge length is zero and the opposite edge length is n 
and l a b € (0, n) then 

s(k) ~ (-l) 2 ^^ 1 / 2 r(fc)~ 1 / 2 sin(U)- 1 / 2 . 

(g) If t exists, but all edge lengths are or n then s(k) = 1. 

Some of these formulas hold for slightly more general sequences of labels, 
see Theorems 14.5. 215.0. 41 below. Most of the results were checked numerically. In 
Figure 3 we show a numerical comparison generated by Maple between the 6j 
symbols and the asymptotic formula Theorem 11.0. II (a) for the sequence of 6j 

symbols I _„ _^ i as j varies from to 108, with r = 200. 2 The phase 



52 54 j 

function <f> was computed by numerical integration, using Schlafli's formula 12.4. II 
(jmjl . The amplitudes ±2nr(k)~ 2 det(cos(/ a b))" are also shown , as well as the 
functions from case (b) of the main theorem governing the degenerate limits . 

An outline of the paper is as follows. Section 2 contains background. The 
conjecture on the geometry of conformal blocks in the semiclassical limit is ex- 
plained in Section 3. Sections 4 to 7 contain a proof of the main theorem. Section 
8 contains a short discussion of the classically forbidden case. The results for the 
hyperbolic case, that is, the case q real and positive, are stated without proof in 
Section 9. Various questions we were not able to resolve are listed in Section 10. 

Acknowledgements The project was suggested by J. Roberts at the end of 
|34j . and we thank him for his encouragement. Discussions with L. Rozansky, F. 
Luo, I. Korepanov, M. Leingang, M. McDuffee and I. Rivin were also helpful. The 
reference |36j and the idea of using recursion to prove the formulas were pointed 
out to us by N. Reshetikhin. Related work on the asymptotics of Turaev-Viro 



2 For the purposes of Maple, dominant weights were labelled by non-negative integers. 
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Figure 3. Quantum 6j symbols versus the asymptotic formulas 



has been carried out by Frohman and Kania-Bartoszynska |15j and J. Murakami, 
unpublished. 

2. Background 

2.1. The quantum group [/^(sb). The quantum group U q (sl2) was introduced 
by Kulish and Reshetikhin in |23j . Let q be a complex number not equal to 0, 1 or 
— 1. For any integer n, the quantum integer [n] is defined by 



Note that [n] — > n as q — > 1, and also 

sin(mr/r) 
M = — :~~7 — r~r~ it g = exp(7r«/r). 
sm(7r/r) 

Let sl2 denote the Lie algebra of 2 x 2 traceless matrices, with basis E, F, H and 
relations [H, E] — E, [H, F] = —F, [E, F] = 2H. The irreducible representations 
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of 5I2 are Vj, j = 0,^,1, ... , where Vj is the space of degree 2j homogeneous 
polynomials in two variables z, w. A basis for Vj is given by Vjj, Vj t j—i, . . . , 
where the weight vector Vj tTn is defined by 

Vj,m — z a w h , a — b = 2m, a + b = 2j. 

The action of the 5(2 in this basis is given by 

Evj,m = (j ~ rn)vj. m+ i, Fvj. m = (j + m)vj tm -i, Hvj. m = mvj. m . 

The quantum group U q (s\.2) is the algebra with generators E, F, K, K^ 1 and rela- 
tions 

K 2 _ K -2 

IE, F] = =— . KE = qEK, KF = q FK. 

q-q 1 

The action of U q (s [2 ) on Vj is defined by 

Evj,m = [j - m]v j>m+i , Fvj !m = [j + m]i?j )m _i, Kv jjm = q m v^ m . 

The coproduct on Ug^sfe) defined by 

E^E^K + R- 1 ®E, F^>F®K + K~ 1 ®F, K ±x i-> K ±l ® K ±l 

makes the category r R,{U q {s{<z)) of finite dimensional U q {si-i) modules a tensor cat- 
egory. TZiUqislz)) has a semisimple subquotient T{U q {s\.2)) called the fusion cate- 
gory for C/ g (s I2) @] Section 3.3]. Any simple object in T{U q {s\.2)) is isomorphic to 
Vj for some j satisfying 



J 



0, 1/2, 1, . . . if q is not a root of unity ; 

0,l/2,l,...,(r-2)/2 if g = exp(f) . 



If g = exp(7ri/r) the branching rule for the product is 

Vj^Vj^QVj 



"3 

3 



where j satisfies the quantum Clebsch-Gordan rules (JIJ. If q is not a root of unity 
the branching rule is that for representations of s^. 



2.2. The 6j symbols. 6j symbols for SI2 were introduced as a tool in atomic 
spectroscopy by Racah .33] 1 and then studied mathematically by Wigner; see |42) . 
6j symbols for U q {sl2) were introduced by Kirillov and Reshetikhin [22], who used 
them to generalize the Jones knot invariant. The material below can be found in 
the book by Carter, Flath, and Saito 10 . If Homy i s i 2 )(Vj, Vj 1 ®Vj 2 ) is non-trivial 

V 

we fix as a generator the map X defined in |1U1 p. 96]. (A small difference is that 
in [101 the map is from Vf j to V? n ®Vf n .) If Hom^ (s[2) (V,-, V jl <g> V h ) is trivial, 
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define Y =0. Let j ao , 1 < a < b < 4 be non-negative half- integers. Define 



I 12 /323 /334 




and 



312 J23 /i 34 



323 J34 312 324. 



312 324 
° 014 



where denotes the identity. By associativity of the tensor product we have a basis 
for Hom^^^V^, Vj 12 ® Vj 23 (g) V^ 34 ) for each way of parenthesizing: 

t} 2 /323 /334 

Hom^ (a , a) (^ 14 , (v jl2 ® y J23 ) ® y j34 ) = 0c 

and 



J13 



Hom^ ( , b) (^ l4 , V jl2 ® (F j23 ® V h J) = < 



324 



where the sum is over ji 3 resp. j'24 satisfying the quantum Clebsch-Gordan in 



equalities for each vertex. The 6j symbol 



12 J23 Jl3 



■334 



\ hi ju hi J q 

' 12 /J23 /334 



is defined by 



E 



Jl2 J23 Jl3 



J34 J 14 .724 I n 

Jl3 




Assume that q is real or |q| = 1, so that the quantum integers are real. One 
can modify the definition slightly so that the symbols have tetrahedral symmetry: 



3l2 323 313 

hi hi hi 



(-1) 



J12+i23+i34+J14 



[2J13 + 1] 



0(123)0(134) 



0(234)0(124) 



jl2 J23 313 

hi 3u hi ) 



where 

q( h \_ ihb + jbc - jac] ! [jab ~ jbc + jac] j [-jab + jbc + jac] j [jafc + jbc + jac + 1] j 
((1 C> ~ (-iy^ + ^[2jab}Wbc]WacY. 
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and the quantum factorial is defined by 
(4) [n]! = [n][n-l]...[l]. 

That is, 

jl2 323 jl3 1 _ / io-(l)cr(2) ia(2)<7(3) 3a(l)a{3) 
.734 jl4 J24 J \ io-(3)cr(4) J CT (1)<t(4) 3a(2)a(4) 

for any cr € S4. Equivalcntly, the 6j symbol is invariant under the exchange of any 
two columns, and any two entries in the first row with the corresponding entries 
in the second. A word on sign conventions: In the case q arbitrary one defines 
the 6j symbols without taking the absolute value of the 9 values. This leads to a 
difference of sign (— l)£a<6 •?<><> from our convention, which was chosen because it 
agrees with the usual definition of classical 6j symbols. For readers familiar with 
spin networks, the 6j symbol is the value of a tetrahedral graph with labels j a b, 
divided by the square root of the absolute value of the product of the thcta graphs 
for the faces. This graph is dual to the edge graph of the geometric tetrahedron, 
that is, edges j a b,jbc,jac are joined at a vertex, for any l<a<6<c<4. 
The Qj symbols satisfy the pentagon or Biedenharn- Elliot relation 



(5) 



E 



(-l)*[2ji B + l] 



J15 



J23 J34 .724 I I J23 .734 ]2A 

hi 312 jl3 J 1 345 325 335 

313 334 314 \ f 312 324 jl4 I I Jl2 J23 Jl3 

345 315 335 J \ 345 jl5 325 J \ ]35 jl5 325 



where z = ji 2 + j\3 + 3i4 + 3i5 + 323 + .724 + 325 + 334 + 335 + 345- This identity is 
obtained by applying the definition of the 6j symbols to Figure 4. 





Figure 4. The pentagon identity 
Racah gave a hypergeometric formula whose generalization to arbitrary q is 
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(6) ( j . 12 j23 Jl3 X = A(123)A(134)A(234)A(124) V(-l) z [z + lll/fz)" 1 

[ 534 5l4 524 J ^ 

where 

( 7 ) f( z ) = \ z - 3i2 - 523 - jia]![« - J13 - J34 - 3i$[z ~ h3 - 534 - 524]! 

[z - jn ~ 324 - ju] ! [j'l2 + 523 + 534 + 514 - 2] ! [il2 + 513 + 534 + 524 - z] ! 

[.723 + jl3 + 514+524 - z]\, 
1/2 

[jab + jbc - iac] ! [jab - jbc + jac] ! [-jab + jbc + jac] ! X 



A(afec) - . ., 

[3ab+3bc+3ac+lV- 

and the sum is over integers z such that the factorials in f(z) are defined. Because 
a large amount of cancellation occurs, JBJ is useful for studying asymptotics only 
in special cases. 

2.3. Asymptotics for q = 1. The asymptotics of the 6j symbols for q = 1 were 
considered by Wigner 05j as follows. Consider a system of particles A,B,C with 
total angular momenta 512,523,534- Given that the total angular momentum of the 
combined AB system is j'13, and the total angular momentum of the combined ABC 
system is J14, the probability of measuring j'24 for the total angular momentum 
of the combined BC system is the square of the 6j symbol. Wigner computed the 
classical probability and arrived at the conjecture 

1 2 - 

512 523 513 I _ 1 



534 5i4 524 J 247tVo1(t) 

where r is the Euclidean tetrahedron with lengths l a b = j a b + 5, if it exists, and the 
approximation is meant to hold after averaging out local oscillations. The Wigner 
conjecture was refined by Ponzano and Regge (32] as follows. Let 9 a b denote the 
exterior dihedral angles of r, and 



ibOab- 

a<b 

Ponzano and Regge conjectured in the semiclassical limit 
( 8 ) J 512 523 513 1 _ , cos(> + 7r/4) 



534 514 524 J v /12ttVo1(t) 

assuming that r exists and is non-degenerate, and also formulas for degener- 
ate tetrahedra with non-degenerate faces and the classically forbidden case. The 
Ponzano-Regge conjectures were further refined by Schulten and Gordon see 
Section 4 below. The mathematical meaning of the semiclassical limit is not entirely 
clear. An obvious interpretation (but perhaps not the only one) is the asymptotics 

of the sequence < /^ 12 /? 23 /^ 13 1 , k — > 00. In this setting the Ponzano-Regge 

I «534 K514 K524 J 

formula (JSJl has recently been proved by J. Roberts [31], see also [TI],[1I]- The con- 
jectures for the other cases follow from the recursion in |3(i| and its mathematical 
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justification by Geronimo, Bruno and Assche [TBI Theorem 3.8]. The remaining 
cases (c)-(f) can be handled using Stirling's formula applied to ©, see also Sec- 
tions 6 and 7 below. 

2.4. Non-Euclidean tetrahedra. Let E n ,S n ,H n denote Euclidean, spherical, 
or hyperbolic n-space respectively. Let a be an n-simplex in E n ,S n or H n . Let 
lab, < a, b < n denote the lengths of the edges, and I the matrix with coefficients 
l a b- The study of which edge lengths occur is called distance geometry in Blumen- 
thal [Hj- Let Go (I) (also denoted Go(cx) denote the Cayley-Menger (n + 2) x (n + 2)- 
matrix (actually a slight modification) written in block-diagonal form 







' 


1 


1 


1 . . 


1 






1 





_I/2 
2 l 01 


2 l 02 ■■ 


_I/2 
2 'On 


(9) 


G (l) = 


1 


2 l 01 





_ 1/2 

2 l 12 ■• 


_I/2 
2'ln 






1 


-1/2 
2 'On 


2'ln 








For a spherical or hyperbolic simplex a we define the Gram matrices G±(l) (also 
denoted G±(a)) by 

G+{l)ab = cos(/ a6 ), G-(l)ab = - cosh(Z ab ) 
and |G±(I)| its determinant. We will need later the following facts: 

Proposition 2.4.1. (a) // an n-simplex a with edge lengths l ao exists in E n , 
then a has volume given by (n! VoI(ct)) 2 = — |Go(/)|. 

(b) An n-simplex a with edge lengths l a b exists in E n if and only if all principal 
minors in Gq(1) containing the W-entry are non-positive. In this case, a 
is non-degenerate if and only if |Go(0 1 < 0. 

(c) An n-simplex with edge lengths l a b exists in S n if and only if G+(/) is 
positive semidefinite. In this case, the simplex is non- degenerate if and 
onlyif\G+(l)\>0. 

(d) An n-simplex with edge lengths l a b exists in H n if and only if all prin- 
cipal minors of G-(l) are non-positive. In this case, the simplex is non- 
degenerate if and only if \G-(l)\ < 0. 

(e) // a tetrahedron r with edge lengths l a b exists, it is unique up to isometry 
(not necessarily orientation-preserving) in E n resp. S n ,H n . If r is non- 
degenerate, there are two equivalence classes of tetrahedra up to orientation- 
preserving isometry, called mirror tetrahedra. 

(f) In the case n = 1 we have 

\G (l)\=-l 2 12 , \G+(l)\ =sin 2 (/ 12 ), |G_(0| =-sinh 2 (Z 12 ). 

(g) In the case n — 2, there are factorizations, 

\G (l)\ = -4s(s - Z 12 )(a - l 23 )(s - Z 13 ) 
\G + (l)\ = 4sin(s) sin(s - h 2 ) sin(s - I23) sin(s - Z13) 
|G_(Z)| = — 4sinh(s) sinh(s — I12) sinh(s — Z23) sinh(s — Z13) 
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where s = (li2 + ^23 + ^i3)/2 is the semiperimeter. (These quantities appear 
in the numerators of the Heron area formulas.) 
(h) A Euclidean, resp. spherical, resp. hyperbolic triangle exists if and only if 
(fT0f> . resp. JTIJl and (JTUJ), resp. hold: 

(10) lac — lab ~l~ Ibci lab — lac Ibc: Ibc — lab ~l~ lac- 



(11) l ab + l bc + l ac < 2n. 

(i) In the case n — 3, a tetrahedron r with edge lengths l a b exists in E 3 resp. 
S 3 ,H 3 if and only if l a b satisfy the triangle inequalities for each face and 

det<3 (0<0, resp.detG + (l) > 0, resp. det G_ {I) < 0. 

(j) In i/ie Zimit o/ small lengths, the determinants of G± approach the deter- 
minant of the Cayley-Menger matrix, that is, 

det G±(el) 
hm = TdetG o (0- 

(k) A spherical, resp. hyperbolic n-simplex with edge lengths I exists if and only 

*/0 < Z < arccos(— 1/n), resp < I. 
(1) The volumes of non-degenerate n-simplices a in E n ,S n or H n satisfy an 

identity due to Schlafii 

(n - l)KrfVol(cr) = Area(F)d(9 F 

where the sum is over (n — 2) -dimensional faces F of the simplex a, 6p is 
the exterior dihedral angle, and k = 0, 1, — 1 resp. is the curvature. 
(m) For non- degenerate tetrahedra r in S 3 , H 3 ,orE 3 , 

2 K dVol(T(s)) =Y^,hbd0ab- 

a<b 

(n) The derivative of dihedral angle with respect to opposite edge length is 

— \\Go(l)\\~ 1 / 2 l c dlab in the Euclidean case 

= < — |G+(Z)| _1 / 2 sin(Z c d) sin(Z a h) in the spherical case 



fd0 at 



\ dl 

^ c ' y — ||G_(Z)|| _1 / 2 sinh(/ C( i) sinh(/ a (,) in the hyperbolic case 

Proof, (a-d) are all discussed in [5]- Note that the non-degenerate cases of (c,d) can 
be treated uniformly using a theorem of Jacobi ^1 p. 303], as was pointed out to 
us by F. Luo. The case of degenerate hyperbolic simplices is slightly problematic, 
and we do not know a reference; it can be proved by induction on n. We will not 
use this case in the paper, (e-j) are all straight-forward, (k) follows from the fact 
that for all edge lengths equal to I we have 

|G+(Z)| = (l-cos(0) n-1 (l + ncos(Z)), = -(cosh(Z) - l) n-1 (n cosh(Z) + 1). 

For (l),(m) see [23 p. 281]. We remark that a formula for the volume of a hyperbolic 
tetrahedron in terms of its edge lengths is given in Murakami-Ushijima |24j . see 
also Mohanty P^J- (n) ||Go(Z)|| refers to the absolute value of the determinant. The 
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Euclidean case is a computation of Wigncr. We prove the claim for the spherical 
case. Let Vi, . . . , V4 be matrices in SU (2) = S 3 , and Eij — ViV^ 1 . We may assume 

pUab 

v b 



v a 



1 
1 



< " 

e- a " b 



Let Dg = diag(exp(«6*/2), exp(— 16/ 2)). Rotation around E a b by angle 9 corre- 
sponds to conjugation by Dg. We compute 

d 



d9 



Ig^TviDgVcD^V- 1 ) = -2det 



V 3 V 3 

' c v d 



V 4 



where V? — 1,2,3,4 are the components of V; as a vector in R 4 , using the 
diffcomorphism SU(2) — > S 3 . Let V be the matrix with coefficients V? . Since 

(V T V) ij =V i -V j = co S (l ij ) = G+(l) ij 

we have det(V) 2 = det(G+(?)). Since V is in block-diagonal form, 



(12) 
(13) 

(14) 



det(V) 



det 



V 1 K 1 

v a y b 

V 2 V? 

"a v b 

1 cos(/ ab ) 
sin(Z afc ) 



= - sin(Z a6 ) sin(Z cd ) 



det 



v c V 3 

v d 4 



v 4 

v c 

ldTr(£ cd )) 



d9 



d9ab 

died 



which proves the spherical part of the lemma. The hyperbolic case is similar, using 
singular values instead of eigenvalues. □ 

The following gives a particularly useful formula for the dihedral angle of a 
spherical tetrahedron. Let t> , ■ • • > v n m F° r anv {Hi ■ ■ ■ ,ik} C {1, . . . , n} we 
denote by <j(i\, i^-, ■ ■ ■ , ife) the n — k simplex spanned by Vj for j ^ i\, . ■ . ,ik- 

Lemma 2.4.2. For a, b e {1, . . . , n} distinct, let 6 a b be the exterior dihedral angle 
between a(a),a(b). n VoI(ct) Vo1(ct(o, b)) = {n - 1) Vo1(ct(o)) Vo1(ct(6)) sin(0 a6 ). 

Proof. Let h a denote the orthogonal projection of v a to onto the orthogonal com- 
plement of the span of o~(a, b) for i = a,b. 

n(n- I)VoI(ct) = ||/i a ||||/i 6 || sin6»Vol(cr(a,6)). 

Using (n— 1) Vol(a(a)) = Vol(a(a, b)) and similarly for a(b) proves the claim. 

□ 

Corollary 2.4.3. Let t be a spherical tetrahedron with vertices v a ,Vb, v c , v d e S 3 . 



(15) I G+ (t) I * sin U = I G+ (r (d) ) 1 2 1 G+ (r (c) ) 1 3 sin ab , 

where G+(r) denotes the Gram matrix with vertices t etc. Furthermore, 



(16) 



2|G + (r(d))|V2|G4r( C ))|V 2 c s^ = ^±^l. 

acos« cc ( 
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Proof. (|15|l is a special case of Lemma 12.4.21 (llfil) is the derivative of (|15J) with 
respect to cos l c d- We start with squaring (fT5|l : 

(17) \G + (T)\sin 2 l ab = \G + (T(d))\\G+(T(c))\sm 2 9 ab . 

Taking the derivative of the left hand side gives 

d\G+(r)\ 2 

-5 ; — sin lab- 

a cos led 

The derivative of the right hand side is 



I G+ (r(d) ) I \G+ (t(c)) |2 sin at> cos fl o6 



d#ah 



dcosZ cd ' 



Substituting l|15|l in 117|) . one obtains 



(18) |G + (r(d))| 1 /2| G+ ( T(c ))|i/2 2c os0 Qb |G + (r)| 1 /2 s i nZat 



d9 ab 



d cos l cd 



We compute 
d0 ab 



d6 a b dl c d 



( sin Z a fc sin I 



cd 



-1 

sin L. 



dcosZ cd dZ cd dcosZ C(i V l G +( r )l 1/2 
Substitute this derivative in the equation l|18|l . we get 

27t\G+(t(cI))\ 1 ^\G + (t(c))\ 1 ^ cos0 Qb sin 2 l ab . 

Thus, the derivative of the whole equation is 

d\G+(T)\ 



d COS led 

as required. 

We remark that 

d|G+(r)| 

d COS Ld 



= 2\G+(r(d))\ 1 / 2 \G + (T(c))\ 1 / 2 cos0 ab , 



1 COs(Z a fc) COs(/ a d) 

cos(Z a6 ) 1 cos(l bd ) 
cos(/ ac ) cos(/ 6c ) cos(Z C d) 



□ 



3. Geometry of conformal blocks 

In this section we explain the geometry of the conjecture 10 ■ Let pi,P2 ,P3,Pi 
be distinct points on P 1 . Given Z a ( a +i) G [0, 7r], a = 1, . . . , 4 let 



ta 











Let .M = ^(^12,^23,^34,^41) denote the moduli space of flat ST/(2)-bundles on 
P 1 — {pi,P2,P3,P4:} with holonomy around p a conjugate to t a . The fundamental 
group of P 1 — {pi ,P2tP3,Pa\ is generated by the loops 7,- around pj , with the single 
relation 71727374 = 1. Hence 

M = {(ei, e 2 , e 3 , e 4 ) G 5[/(2) 4 , eie2e 3 e4 = /, e Q ~ i„, a = 1, 2, 3, 4}/S77(2). 
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Under the diffeomorphism SU(2) — ► S 3 , the set of matrices conjugate to t a is 
mapped to the set of points in S 3 at distance l a i a +i) from (1,0,0). The 4-tuple 

v 1 = I, v 2 = ex, w 3 = exe 2 , w 4 = eie 2 e 3 

maps to the vertices of a closed 4-gon in S 3 with edge lengths l a ( a +i)- From now 
on, we assume l a t a +i) < 7r - Then .M = -M(7 12 , 123, I34, hi) is the moduli space of 
closed 4-gons in S 3 with edge lengths h( a +i), where we identify two 4-gons if they 
are related by an isometry of S 3 . Ai is non-empty if and only if the inequalities 

hi < I12 + hs + hi, hi + I23 + hi <2n + hi 

and their cyclic permutations hold; this can be seen either directly or from 0,13 
by computing Gromov-Witten invariants of P 1 . If one of these equalities holds 
with equality, then the 4-gon is contained in a geodesic on S 2 and M. consists of 
a single point. If one of the lengths is zero then M. is single point. Otherwise, M. 
is homeomorphic to S 2 . If in addition ±Zi 2 ± I23 ± hi i hi is not a multiple of 2ir 
for all choices of sign then M is non-singular and symplectomorphic to S 2 . 

There are two canonical coordinate systems on Ai, described as follows. Let 
A13, resp. A24, denote the distance from the vertex v\ to V3, resp. v 2 to v^. In 
the language of flat bundles, the functions A13, A24 measure the eigenvalues of the 
holonomy of the connection around a circle separating the four points into two 
groups of two. Let 6*13, resp. #24, denote the corresponding dihedral angle. Let Ai* 
denote the subset of 4-gons that are not contained in a geodesic, that is, the subset 
of irreducible bundles. Let 

M* 13 CM*, M* 24 cM* 

denote the subset on which < A13 < n and the faces containing the edge V1V3 
are non-degenerate, resp. < A24 < tt and the faces containing v 2 v 4 are non- 
degenerate. A special case of a result of Goldman states that Ai 3 ,0i 3 , re- 
spectively A 2 4, #24; are smooth on Mi 3 , resp. A^24 and form a set of action- angle 
coordinates. That is, 

u\m$ 3 = d#i3 A dAi 3 , u\ M ^ = d6>24 A dA 2 4- 

If I12 — I23 and ^34 = Z41 then A^ L3 1 (0) is a subset of M. homeomorphic to S 1 , 
meeting the complement of M* in the two configurations where all vertices are 
colinear. Let If 1 ™ denote the minimum value of A13. If Z™ 11 ^ 0, the set defined by 
A 13 = lfi n is a point. If ^ 3 in = 0, the set A 13 = ^ 3 in is a circle. 

By a theorem of Mehta-Seshadri [21] and Furuta-Steer [Tfi] . M has the struc- 
ture of a normal projective variety. By a theorem of Pauly M. has a positive 
line bundle C whose sections H°(Ai, C r ~ 2 ) may be identified with the space of 
genus zero WZW conformal blocks at level r — 2 with weights //x, M2, M3j M4- Let 
Co, 4 denote the moduli space of stable 4-pointed genus zero curves; the open sub- 
set Co, 4 of 4-pointed smooth curves is an open subset identified with C — {0, 1} by 
the map (P 1 ; 0, 1, 00, z) > z. Let X13, resp. X24, denote the point in Co, 4 repre- 
sented by a nodal curve with two components and points p± , p 2 on one component 



16 



YUKA U. TAYLOR AND CHRISTOPHER T. WOODWARD 



and pz,Pi on the other, resp. pi,Pi on one component and P2,P3 on the other. 
Let V denote the bundle with fiber V(ri tPliP2t p 3tP4 ) = iJ°(£ r_2 (P 1 ,pi,p2,P3,P4))- 




According to the Khizhnik-Zamolodchikov, see [HHjj|S]j the bundle V has a projec- 
tively flat connection V. Furthermore, the fiber over X13, resp. X24 has a canonical 
projective basis {sis(ji3)} resp. {324(^24)}- The 6j symbols are the entries in the 
change of basis matrix for the bases {si3(ji3)}, {S24O24)} after parallel transport 
to a common point in 6*0,4. (However, we only know how to give a canonical 
normalization to the basis elements using the quantum group definition.) Define 

A» - (£5) , A, - A- . 

A13, A 2 4 are Bohr-Sommerfeld Lagrangians. That is, the restriction of C r ~ 2 to A a b 
has a flat section s^O'ab)- Extend the flat section to a section s^(j a b) defined 
in a neighborhood by parallel transport in the perpendicular directions, followed 
by multiplication by the Gaussian function exp(— 4n(r — 2)d 2 b ) where d a b is the 
geodesic distance to A a {,. Soo(jab) is not in general holomorphic. We conjecture 
that as k — > 00, after replacing (r — 2) with k(r — 2) one has 

(19) s ktab (kjab) ~ s^ ab {kj ab ). 

Part of the problem is to give a precise sense in which this asymptotics holds. 

We briefly discuss a possible approach to proving lfiT?|) . The vector field ^J^- 
has (where defined) a canonical lift to a connection-preserving vector field on C 
However, ^J^- does not preserve the holomorphic structure, and so even if globally 
defined does not define a circle action on H°(£ r ~ 2 ). However, we conjecture that 
in the limit (P 1 ,pi,p2,P3,P4) — + A13 in Co, 4, the complex structure is invariant; 
see Daskalopolous-Wentworth ^2]- This would imply that the Verlinde sections 
S13C713) are asymptotically eigenvectors for the circle action with weight j'13. I|19|) 
would follow from standard arguments discussed in We also remark that the 
leading order part of V as discussed in ^3, (4.23)]) is a first order equation and 
should not change the large k asymptotics. 

The determinant line bundle C is conjectured to admit a Hermitian metric 
so that the induced Hermitian metric on V is compatible with the projectively 
flat connection [3], see also 031 jEJ- It would follow that the 6j symbol equals the 
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Hcrmitian pairing of the sections (S13C713), S24(,?24))- Since Sk,ab(kj a b) are asymp- 
totically Gaussian concentrating to A a h, one naively expects the pairing to have 
asymptotics 

(2Q) J2 {k[r - 2 ))' 1 (^M3(fcji 3 )(r), Sfc , 2 4(fcj 24 )(T)) e ± "/ 4 _ 

reAi 3 nA 2 4 \J ( kuJ )r(g§^, g§^) 

We denote the two points of intersection t + ,t_. The pairings w T± (jj|— , are 
bv EXTl tatl 

( d d \ dl 13 _ |G+(0l 1/2 



\B9i3 862aJ d024 sin(/z c d) sin(^ fc) 
Moving along the manifold A13 from r + to t_ and then A24 from t_ to t + produces 
a path 7 shown in bold in the Maslov diagram in Figure 2, case (a). Let Area(7) 
denote the symplectic area enclosed by 7. Then 

(Sl3 (jl3 ) , S 2 4 C?24 ) ) (t+ ) = CXp (i Aie'd,(j) ) (si 3 (j'13 , S 2 4 ( J24 ) ) (t_ ) . 

(I2L)[1 becomes 

\ v, cos(iArea(7) + f) 

l(Sfc,13(« : Jl3),Sfc,24(Kj24)| ~ , =■ 

(fc(r-2))3/ 2v /c(^,^) 

We claim that 

Area(7) = 2 Z afc 6> afc + 7r(/ 12 + Z 2 3 + ^34 + hi)- 

a<b 

As we vary Z13 and £24, 

dArea( 7 ) = 2(0i 3 dZ 13 + 8 24 dl 2i ) = 2{d^0 ab l ai - 2 Vol(r)) 

a<b 

by Schlafli's identity 12. 4. II ifmjl . Hence 

Area(7) = 2 ^ 6» ab / afc + c(Zi 2) I23, hi, hi) 

a<b 

where c(l\2, I23, I34, In) is a constant depending on £12, £23, £34, hi- For Z13, £24 cor- 
responding to a degenerate tetrahedron with four exterior edges, we have 

#a(a+l)=7T, a =1,2, 3, 4, 13 = 6» 24 = 0. 

This implies 

c(£l 2) ^23, ^34, ^14) = 2tt(Zi2 + £23 + ^34 + £14)- 

Putting everything together, this gives in the transversal case the naive prediction 

27rcos (*^(£ W«4 - 2 Vol(r)) + f ) 
|( S M3(^i3), S M4(fc j2 4))| ( fc(r -2))3/ 2 det(cos(U)) 1 / 4 ■ 

Including half-forms in the quantization scheme should produce shifts giving Theo- 
rem ll.U.fl fa'). (Readers familiar with the representation theory of affine Lie algebras 
will recognize the numbers ^,2 in the formulas for l a \, and r(k) as the half-sum 
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of positive weights, resp. dual Coxeter number of s( 2 .) Unfortunately we lack the 
techniques to carry this out. In the follow section we introduce a recursion which 
reduces the formula to the q — 1 case; the recursion is needed anyway to handle 
the first degenerate case, so even if a geometric proof could be found for case (a) 
the proof of the Theorem II. 0.11 we give is presumably shorter. 



4. A NON-DEGENERATE TETRAHEDRON 

First we show that both sides satisfy a second order difference equation as 
one label is varied; the same strategy was followed by Mizoguchi-Tada [27| without 
knowledge of the denominator in the asymptotic formula. This implies that each 
side is a linear combination of the two linearly independent solutions. To prove 
that the coefficients are equal we show that in the limit as the labels jab a re made 
small in relation to the level r — 2 our equation reduces to the Ponzano-Regge 
formula JSJ). 

4.1. The Schulten-Gordon recursion. 
Proposition 4.1.1. 

(2i) p{m){f f f A } + [2j 23 ]Q(; 2 3 + 1) { f J23 7 + 1 f 

I J34 Jl4 J24 I I J34 314 324 



[2j 23 + 2]Q(j 23 ) 



jl2 J23 - 1 J13 
J34 jli J24 



where 

P(j23) = [2j 23 ] [2j 23 + 1] [2j 23 + 2] \j U + J 13 + J34 + 1] [,?13 + 334 ~ hi] 

- [2J 23 ] [jl2 + jl3 - j 23 ] [jl2 + J23 ~ jl3 + 1] [?24 + J34 ~ J23] [hi + J23 - 334 + 1] 

- [2j23 + 2] [j 12 + J23 + 313 + 1] [jl3 + J23 ~ J12] [j24 + J34 + J23 + 1] L?34 + 323 ~ 32a] 

and 

QU23) = ([jl2 + J23 - 313] [jl2 + 313 - 323 + 1] [324 + 334 - 323 + 1] L?24 + J34 ~ j 23 ] 
[jl2 + j23 + 313 + 1] [jl3 + 323 ~ jl 2 ] [j24 + j 3 4 + 323 + 1] [j34 + 323 ~ J2i\ ) 5 ■ 

The proof is a combination of pentagon identities JSJ, as in [36) . and is omitted. 

In the limit r — > oo the Schulten-Gordon equation (|23|l simplifies dramati- 
cally. Let A ab denote the discrete Laplacian 

(&abf)(jab) = ftiab + 1) - 2f(jab) + ftiab ~ !)• 

Let t denote the tetrahedron with lengths l a \, = (j a {, + ^)/r, if it exists, and Q a \, 
the exterior dihedral angle around the edge e a fc. Set 

i 

ff j s_( \ G +Q)\* ) 2 / jl2 J23 jl3 
/Uabj " I sin(M J \ 334 314 324 
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Let Ts C K 6 denote the subset of edge lengths of possibly degenerate tetra- 
hedra, with all faces non-degenerate. Ts is defined by the inequalities 

(22) \G+(T(d)) > 0|,1 < a < b < c< 4 

|G+(r)| >0. 

Theorem 4.1.2. Let K be a compact subset of Ts and 1 < a < b < 4. There 
exists a constant C and a second order difference operator e(r) with coefficients 
bounded by Cr~ 2 for r sufficiently large, such that 

(23) (A afc + 2 - 2cos0 ai , + e(r))/0'o&) = 
on the relative interior of K . 

Remark 4.1.3. One may extend this recursion to the subset of R 6 of 6-tuples l a t 
satisfying only l|22[) by defining (9 ab to be the imaginary number given by (|16|) . 

Proof. Let cd be the opposite edge of ab. A somewhat long computation involving 
standard trigonometric identities shows that 



(24) 



< Cr- 2 . 



Substituting this into the recursion Q21|l and using the factorization l2.4.Tl Q gives 

(25) lMi ab )^p\l j r 3ab f ) + - A Mi ab --y 

y ' 4 'dcos(Z cd ) ^ Jbd Jed Jad J 4 r 

2ir 2tt i 

[\G+(l a c, lab + — , lbc)\\G + (lbd, lab + — , lad) \ \ G + (l ac , l a b, he) \ \G+ (hd, lab, lad)\]^ 

■<:- jab 1 1 jbc Uisin(u + -)-| jac 1 •"" 

Jbd Jed Jad J 4 r { Jbd Jed Jad 

2tt 2tt i 

\\G+{l a c,lab ,lbc)\\G + (lbd,lab , hd)\\G + (l ac , hb, hc)\\G + (hd, hb, hd)\] 7 = 0(r ) 

r r 

on the relative interior of K. (Here G+{l a c,lab,hc) etc. refer to Gram matrices of 
the faces.) We divide by 

i 

2tt 2n \ 2 i 

sin(; ab ) sin(; ab ) sin(; ab H ) [\G+(l ac , l ab , l b c)G + (l a b, hd, Ld)\Y ■ 

r r J 

Using the equality 

sin(U)(sin(U - 2vr/r) sin(/ Qb ) sin(U + 2njr))^ = sin(U) _1/2 (l + 0{r~ 2 )) 
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and similar identities one obtains the vanishing of 
(26) 

~sm(l ab )-^\G + (l ac J ab J bc )\\G+(l ab ,l bd ,l ad r~- dlG+(l)l I Jac Jab Jbc 



dC0S(7 cd ) !_ jbd jcd jad 

+sm(l ab +— )- 1 mG+(l ac> l ab +— ,l bc )\\G + (l bd) l ab +— Jac Jab + 1 Jbc 

r r r y ]M Jcd Jad 

+sin(U-— rHG+ilacU- — ,hc)G + (l bd ,l ab - — { JaC 3ab ~ 1 ]bc 

r r r [ Jhd Jcd Jad 

up to terms of order r~ 2 . Substituting (|16(l and i|15|) into the above gives the 
result. □ 

Remark 4.1.4. An examination of the proof shows that the constant C depends 
linearly on a lower bound for (\G + (l ac ,l ab ,l bc )\\G + (l ab ,l bd ,l ad )\)i . 

4.2. Second-order recursion for the asymptotic formula. We now show that 
the right-hand side of J2J| is also an approximate solution to Q23p. Set 



Jl2 J23 Jl3 



27rcos( 



(27) 

' J34 314 J2A J <j=exp(7ri/r) r*|G+(/)|* 

Let A^ b be the second-order discrete Laplacian at step size 27r/?', 

(A^/)(I) = f(l - 2n/r) - 2/(0 + f(l + 2n/r). 

Let Tg" denote the subset of R 6 defined by the inequalities (|22Jl and |G+(t)| > 0, 
that is, the set of edge lengths of non-degenerate tetrahedra. 

Theorem 4.2.1. Let K be a compact subset of Tg . There exists a constant C 
and a second-order difference operator e(r) with coefficients bounded by Cr~ 2 for 
r sufficiently large such that 

(K b + 2 - 2008(4*) + e(r)) (^^) ' { f f f V 

y Sm6 ab J { J34 314 J24 ) q=cMm /r) 

vanishes on the relative interior of K. 
Proof. Let 8 = 9 ab . By Taylor's theorem 

(28) <p(l ab ±— )=</>(lab)±-5r- — i.U) + ^^2 — (U + e) 



6 



for some e G [— ^r,^ 2 -]- Because is compact, ^ ^ 3 is bounded on an open 
neighborhood of K. Bv l2.4.ll (jn*T(l . 



r 



d 2 d> r 89 



dU 2tt ' dl ab 2 2-Kdla 
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Thus 

(29) cos(0(U ± 27r/r)) = cos(0(U) ± 0) - -^-sin(0(U) ± 0) + 0(r" 2 ). 

r ol ab 

Similarly 

S in-V 2 (e(l ab ± 27r/r)) = sin- 1 / 2 ^^)) T - sin- 3 / 2 (0(/ o6 )) cos(9)-^- + 0(r~ 2 ). 

r ol ab 

Expanding cos((/>(Z a b)±0), sin(0(i a h)±0) one finds that the 0(1) and 0(r _1 ) terms 
cancel, which completes the proof. □ 

Remark 4.2.2. An examination of the proof shows that the constant C depends 
linearly on an upper bound for ai °| 



dl a 

4.3. Euclidean limit of the 6j symbol. Next we investigate the limit of both 
sides of J3J) as the level is taken to infinity much faster than the labels. 

Proposition 4.3.1. Let j a f, S [0, (r — 2)/2] and 5 a b be half-integers, for 1 < a < 
b < 4, and p > 3 an integer. Let j a b(k) = kf a b + S ao . Then 

ju(k) J23(fc) jl3(k) \ f ju(k) jz3(k) j 13 (k) 



J34(fc) Jli(k) J24(fc) J g=cxp(Wr(fcp)) I J34(fc) Jl4(fc) J24(fe) Jq=1 

The proof relies on the 

Lemma 4.3.2. There exist C, c > swc/i iftcrf ifr/n < c then | In ([n]!/n!) | < 
Cn 3 r- 2 . 

Proof. Applying the logarithm to gives 



ln[n]! = ^l n (sin(^))-ln( £ 

3=1 



sm I — 

r I I \ \r 



Note that | sin(x) — x\ < %- and ln(sin(x)) = ln(x) + (sin(x) — a;)/?/ for some 
y S [x,sin(x)]. For x sufficiently small, 

x 3 

sm(x) €[x- —,x], y>x/2 
6 

which implies 

x 3 2 x 2 

| ln(sin(x)) - ln(x)| < — - = — . 

X 6 

Hence 

| ln[n]l - lnn!| < £ -(( + {-?) < ^((—f + (-?) < 



3 r r 3 r r 3r 2 

3=1 



□ 



Proof of \4-3-l\ By Lemma 14.3.21 and Racah's formula ©, since the number of 
terms in the sum grows linearly with k. □ 
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4.4. Euclidean limit of the asymptotic formula. In this section we compare 
the quantities 

f 312 323 jis 1 °° . = 2ttcos(^+ f) 

\ 334 ju 324 J 9=exp(wi/r) ' ri\G+(l)\i 

and 

f 512 323 313 \°° 27TCOS(0O + f) 

I 534 il4 524 ] q=l ■ \\G (3 + m i/4 

in the case that the corresponding tetrahedra are non-degenerate. We will prove 

Proposition 4.4.1. Let j ab e [0, (r — 2)/2],l < a < 6 < 4 be such that a 
non-degenerate Euclidean tetrahedron with edge lengths j a b exists. Let 5 a b be half- 
integers, and p > 3/2. Let j a b(k) — kj a b + Sab- 

( 512(A) j 23 (A) 513(A) 1 °° f 512(A) j 23 (A) 513(A) 1 °° 

\ 534(A) 514(A) j M (k) j 0=exp(7ri/r(feP)) 1 534(A) 514(A) 524(A) j q=1 ' 

Proof. ByEXll§ r(fcP) 3 / 2 |G+(0| - ||G (j + |)||. Since the edge lengths are 
0(fc 1_p ), the volume Vol(r) is 0{k 3 ~ 3 P). Let 0f 6 denote the dihedral angles of the 
Euclidean tetrahedron with edge lengths l a b, if it exists. It follows from Corollary 
EUthat cos(9 ab ) = cos(0f b )(l + O(fc 2 - 2 P)). Since 9 ab 9% as k -► 00, a6 -0f 6 < 
Ck 2 - 2 P. Hence _ 

which implies that 

0(fc)-0 o (A) <Ck 3 - 2p . 
For p > 3/2, this approaches zero as fc — * 00. □ 

4.5. All edge lengths equal. The second order difference equation corresponds 
to a system of first-order difference equations in the standard way. In fact, if v is the 
vector formed from all 2 6 sets of labels of the form j a b+S a b, where d a b = or 1, then 
the vector of the corresponding 2 6 6j symbols satisfies a first order recursion in all 
six variables jab . We call the first-order systems corresponding to (|21|l . (|23|l the first- 
order Schulten-Gordon difference equation and first-order asymptotic difference 
equation, respectively. By Theorems 14.1.21 and 14.2.11 to show Theorem ILL). II (a) 
it suffices to prove the claim for all components of v and a single set of values of 
jab- We do so for the case that all edge lengths are equal. Bv l2.4.1l (|kjl. a spherical 
tetrahedron with edge lengths equal to I exists for < I < arccos(— 1/3). Let 
h = h(k) be a function of k so that h > k p , for some p > 3. Choose r >> so 
that there exists j S (0, (r — 2)/4) n Z/2. For sufficiently large k, there exists an 
approximation Jk{t) to the linear path lt,t £ [k/h(k), 1] such that 

(a) Yfe(t) + S a b/r(h(k)) is the set of edge lengths of a non-degenerate spherical 
tetrahedron, for any choice of <5 ao ; 

(b) jk is a lattice path for Z/2r(h(k)), that is, obtained by concatenating paths 
between lattice points whose difference is a standard basis vector divided 
by 2r(fc(*)); 
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(c) lk (k/h{k)) ab = {kj + \)/h{k) and 7it (l) o5 = (h(k)j + \)/h{k); 

(d) the number of lattice points in 7^ is 0(h(k)). 

Let l(k, to), to = 0, . . . , n(k) be the lattice points occurring in 7. Let ^4 m be the 
operator in the 771-th term of the first-order Schulten-Gordon difference equation, 
and the operator in the m-th term in the asymptotic difference equation. On 
an open cone containing I, we have an estimate 

(\G + (l ac ,l ab ,l bc )\\G + (l ab ,l bd ,lad)\)* > C\\l\\ 2 . 

By Remark 14.1.41 Az 1 A™ = I + e(m) where e(m) has coefficients bounded by 

(30) Ch- 2 \\l(k,m)\\- 2 . 

Let 

7j(fc, to) = A(k, m)A(k, m - 1) . . . A(k, l)v(k, 0) 

denote the solution to the Schulten-Gordon difference equation, and 

v°°(m) =A oo (m)A oo (m-l)...A oc (l)v{0) 

the solution to the asymptotic equation. Then 

v°°(m) = A(m)(I + e(m))A(m-l)(I + e(m-l))...A(l)(I + e(l))v{0) 

= (/+ e'(m))(J + c (to - 1)) . . . (J + e(l))«(m) 

where e'(l), . . . ,e'(m) are operators with coefficients also bounded by H3U|) . The 
norm of the product of error factors satisfies the estimate 




e'(m)) 



n(fc) 

< C* ^ ||W(fc,m)|| 

m— 1 



which goes to zero as k — > 00. It follows that the components of tj°°(to) are 
asymptotic to those of v(m) as k — > 00. (The convergence of approximate solutions 
of second order difference equations to actual solutions is discussed in much greater 
generality in [T%|.') 

A similar estimate holds for the right-hand side of On the cone generated 
by a small neighborhood of l® b we have a bound 

Olab 

Using Remark 14.2.21 the error factor approaches 1 as k — > 00, that is, the asymp- 
totic 6j symbols approach an exact solution of the asymptotic difference equation. 

To show that Vk{n{k)) is asymptotic to v^(n(k)) it remains to show that 
each component of Ufc(O) is asymptotic to the corresponding component of tj^°(0). 
By Propositions 14. 3. 11 14.4. II it suffices to show that for any 5 ab € {0, 1}, 

kj + 5i2 kj + S 23 kj + £13 1 f kj + S 12 kj + 823 kj + 813 



kj + S 3i kj + Su kj + S 2i J l kj + S 3i kj + 5u kj + 5 2 4 , x 

For 5 ab = this is Roberts theorem |34j . To handle the remaining cases, we need 
a generalization: 
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Theorem 4.5.1. Let ji2(k), . . . ,j 3 i(k) € Z/2 be functions of k such that l a b(k) = 
27r(ji2(fc) + \)/k converge as k — » oo to l a b G (0, tt) the edge lengths of a non- 
degenerate Euclidean tetrahedron. 



oc 



iia(fc) J23(k) ji 3 (k) j ^ j j 12 (k) j 23 (k) j 13 (k) 
j 3 i{k) ju(k) ]2i{k) ) q=l \ j 3i {k) j u (k) j 24 (k) J g=1 

Proof. It seems that this follows by the same method used by Roberts [31] . Alter- 
natively, the 6j symbols for q — 1 as a function of any j ab satisfy the Schulten- 
Gordon second-order difference equation up to an error term depending linearly 
on an upper bound for | |G (/ a fc, Lc, hc)\ \G (l a b, lad, hd)\ |~ 1/4 . Let Tj£ be the sub- 
set of K 6 consisting of edge lengths of non-degenerate tetrahedra, and K C Tg 
a compact subset. On kK, the 6j symbols solve the difference equation up to an 
error operator of order 0(k~ 2 ), and so equal an exact solution up to 0(k~ 1 ). Any 
exact solution is a linear combination 

2irr-V 2 \\G (l)\\- 1/4 ( Cl cos(0 o + tt/4) + c 2 sin(0 o + tt/4)); 

We wish to show that c\ = 1, C2 = 0. Let ab — 12 and suppose that j a b(k) = kj a b 

for ab ^ 12. Roberts' theorem applied to { f j23 k J ls \ gives 

I k ] 3 A KJ14 KJ24 J 

(31) ci cos(0o + tt/4) + c 2 sin(0 o + tt/4) = cos(0 o + tt/4). 
Roberts' theorem applied to ( fc C?i2 + *) Jtfra Jtfia 1 iyeg 

[ fe?34 KJ14 KJ24 J 

(32) ci cos(0 o + 6*i2 + tt/4) + c 2 sin(</» + 0ia + tt/4) = cos(</- + On + 7r/4). 

Since #12 € (0, w), the equations J2U,JS2l are linearly independent, hence c\ = 1 
and C2 = 0. Applying the same argument to the variation in j 23 ,ji 3 etc. proves 
the theorem. □ 

This completes the proof of Theorem 11.0. II (a) . The argument used for Theorem 
14.5. II applies to the 6j symbols for q = exp(-7ri/r) and gives the following general- 
ization of Theorem ll.0.11 (a). 

Theorem 4.5.2. Let r > 2, j\ 2 {k), . . . , j 34 (k) G [0, (r - 2)/2] n Z/2 be functions 
of k such that l a b(k) — 2Tr(j 12 (k) + h)/r(k) converge as k — > oo to l ab £ (0, ir) the 
edge lengths of a non-degenerate spherical tetrahedron. 

jn{k) j 23 {k) j 13 {k) \ ( j 12 (k) j 23 (k) j 13 (k) 



Mk) ju(k) Mk) J , =oxp(ffi/r(fc)) I J34(fc) ju(fc) J24(*) j , =oxp(Wr(fc)) 
5. A DEGENERATE TETRAHEDRON WITH NON-DEGENERATE FACES 

We apply the difference equation of the previous section to a sequence of 
labels in which the corresponding tetrahedra degenerates to a degenerate tetrahe- 
dron with non-degenerate faces. The two possible limiting cases are shown below in 
Figure 6. Let l° ab , 0^ b denote the lengths, resp. angles of the degeneration. Following 
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Figure 6. Two limiting cases 

Schulten and Gordon we define 
_ r{k) ^ „ 



9° 
'ah 



a<b 



0, iffl Qb <^/2; 
7r, if 9 ab > ir/2 . 



In the two limiting cases shown above, we have 

3 

00 = 7r0'l2 + J23 + J34 + J41 + 2), resp. 7r(ji2 + j 2 3 + jn + 

Hence 0o is an integer, resp. half integer, times it. By Schafli's formula 12.4. II (fTTTfl 



(33) 



d(^-0o) 



r(k) 
~2^T 



5^(#afc - ^Qb)d^a6 



which implies — 4>q < 0, resp. — 0o > for tetrahedra near the first, resp. 
second, limiting case. 

We apply the LG-WKB method to arrive at a conjectural formula. In order 
for / to be a solution to the recursion relation l|23|l. fa, fa must satisfy 

(A ab + 2-2 cos(0 ab - 6° ab ))f(j ab ) = 0. 

According to ansatz described by Schulten and Gordon, in the semiclassical limit 
the solution must solve the differential equation 



(34) 



d 2 



djal 



- {Oab - Bab)' 



Let F(x) denote a solution to the Airy equation F"(x) = xF(x). Define 



f{jab) 



sHOab ~ 9° ab ) 



-1/2 



A(j ab )F(n(j ab )). 



Inserting this into 1)34(1 gives 

^ + n' 2 n + (8 ab - e° ab ) 2 ^j F(si(j ab )) + \2^n' + F'(si(j ab )) = o. 

Assuming that A" /A « 0, the solution is 



tt(jab) 



{Oab - 6ab)djab 



2/3 



'Mjab)-0°(jab)) 2/3 



A(j ab ) = C\<fa ab )-4>°Uab)\ 1/ *8£ /2 . 
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We need to choose F, f to agree with the solution in the transversal case. Introduce 
the regular and irregular Airy functions 



If 00 1 f°° 

Ai(x) = - / cos(i 3 /3+xt)dt, Bi(x) = - / (exp(-i 3 /3+rf)+sin 

7T Jo 7T J 



(t 3 /3+xt))dt. 



For large x one has the asymptotics 



(35) Ai(-x) ~ 7r" 1/2 a; _1/4 cos(£ - Bi(-a;) ~ tt- 1 / 2 ^ 1 /^- sin(£ - J)) 
where £ = |a; 3 / 2 . It follows that the function 

( 36 ) 2tt 3 / 2 Z 1 / 4 [cqs(0q) Ai(-Z) - sin(0 o ) Bi(-Z) - O < 
r 3 / 2 sin(0 ah ) 1/2 \cos(<fo)Bi(-Z) - sin(0o) Ai(-Z) - O < 

where 

^ = (^-0o|) 2/3 

satisfies the asymptotic recursion and matches up with the non-degenerate solu- 
tion. 

Lemma 5.0.3. In the limit 9 ab — * 0° ab , 

6 - 6 3r 



|G + (Z)I 3/2 2ttA 1 A 2 A 3 A a 
Proof. Using itTBT) . at l ab — l ab we have 

g |g+(OI 3/2 3 m i_»_, r m.i/2 



|G + (Z)| 



1/2 A a A fc sin(6> ab ) d A e A rf sin(6> erf ) 



AiA 2 A 3 A 4 sin(7 afc ) dl ab sin(7 cd ) 

S\G+(l)\^ 2 A a A b sm(0 ab ) A c A d cos{e cd ) d9 cd 
AiA 2 A 3 A 4 sin(7 ab ) sin(/ cd ) dl ab 

3|G + (0r 



|1/2 sin(6> afc ) cos(9 cd ) sm(l ab ) sin(Z cd ) 



sm(l ab )sm(l cd ) |G+(0| 1/2 
= 3cos(6» C(i )sin(6U). 

The Lemma now follows from L'Hopital's rule. □ 
The values of the Airy functions at zero are 

Ai(o) = (3 2 / 3 r(^))-\ Bi(o) = (3*r(|))- x . 
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Therefore, the limit of the 6j symbol as j a b — * 3%, is (— l)J 12 +J 2 3+J34+j4i reS p_ 
(_l)3ia+j23+jia times 

J^w^f^ = A r»/>iG + (or — Al(0) 



= 2rr 3 / 2 (3/2) 1 / 6 — -4— — 3- 2 / 3 r ± 

= r- 4 /3 2 2/3 3 -l/3 7r 4/3 (AaA() ^ ) -l/6 r /0 1 

which is (b) in Theorem ll.0.11 

The LG-WKB method for finite difference equations near a turning point is 
studied by Geronimo, Bruno and Assche in 18 , Theorem 3.8]. The assumptions 
18, p. 106] that guarantee that an approximate solution converges to the solution 
given by the ansatz translate to the conditions 

(37) (JL^j (9 ab - 9° ab ) 2 £ C°((0, rr)), i = 0, 1, 2, 3; 



( 38 ) (w-) l {0 j b J}o )2 eg°((0,7r)), i = 0,l,2. 

To see that these conditions holds, note that by (|16[) cos 2 (6* a ;> — #° b ) is a smooth 
function of l a b- Hence so is (0 a b — #2(,) 2 for a b sufficiently close to 0, 7r, which proves 
(|3"?|) . (|3*5|) follows since {9 a b — S Q ab ) 2 has a simple zero at l a b — l^b- Condition (iii) 
on |181 p. 106] is achieved by choosing the sign of the parameter t appropriately, 
for t = ±(l ab — lab) sufficiently small. 

Unfortunately we do not know whether there are any configurations of the 
type described in Theorem II .(J. II (b), that is, degenerate tetrahedra with non- 
degenerate faces and all edge lengths commensurate with tt. However, the following 
generalization of Theorem 1 1 . . 1 1 ( h) is non- vacuous: 

Theorem 5.0.4. Let j a b(k) £ [0, (r — 2)/2] n Z/2 be a sequence of labels such 
that l a b(k) = (jab(k) + h)/f(k) converge to the edge lengths l a b of a degenerate 
tetrahedron with non- degenerate faces and l ab (k) — l a b = 0(k~ 1 ). The sequence 

^} 2 )fl ^? Jfl ^ ,,\ 1 is asymptotic to the expression in Theorem\lA)A\ (b) 
j 34 (fc) ji 4 (k) j 24 (fc) J y y l( J 

as k — > oo . 

Proof. By the discussion above, the sequence has asymptotics given by (|3*6^l . By 
Schlafli's formula and Taylor's theorem, 0- 0° < C(9 ab - 9° ab ) = Oik- 1 / 2 ). Hence 
the argument in the Airy function in (|36|l goes to zero, which completes the proof. 

□ 
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6. One face degenerate 

If one face, say 123, is degenerate then there are two possibilities: either one 
lengths, say j 13 , is the sum of the other two, j 12 + j 23 , or j 12 + j 23 + 323 = r - 2. 
In the first case the Racah sum ^ has a single term 

3l2 323 312 +323 \ / [2jl2]![2j 2 3]![jl2 +.?23 + J34 +.714 + 1]! 



(39) 

^_iyi2+j23+j34+jl4 



J34 jl4 J24 J U 2 Jl2 + 2 J23 + l]![-jl2 - J23 + J34 + J14]! 

[jl2 + J23 + J34 ^ jl4]![jl2 + J23 ~ j 3 4 + jljjj 
[.]23 + J34 - J24]![j23 - J34 + J24]![j23 + 334 + J24 + 1]! 

[j34 + J24 - J23]![jl4 + J24 - J12]! 



\jl2 + 324 + jli + 1]![.?12 - jli + J24]![jl2 + jli ~ J24] ! . 

The asymptotics of the 6j symbol are therefore determined by the asymptotics of 
quantum factorials which have been investigated by Moak [21], see also |HJ]. The 
following can be derived from [2EI (2.12-16)]. 

Proposition 6.0.5. (q-Stirling) Set q — exp(iri / kr) . Then as k — > oo, 
(40) [kn]\ ~ V2n[kn] kn+ il(mr/r)- kn 

where 



I(x) :— sin(a;) exp ^— x 1 j 



ln(sin(y))rfy I , lim I{x) 



We conjecture that the factors involving I(mr/r) cancel; in the case q = 1 
these factors are replaced by factors of e which do cancel. Assuming this, one 
obtains 

r ii \ 1/2 

sin(Zi 2 ) ^ sin(Z 23 ) 2 1(li3, hi, hi) 



(41) s(k) ~ 2- 1 ^r(k)-^n 
times (-l) fe 0'i3+J34+ji4) w h ere 

li^ab-i Iba lac) 



sin(Zi 2 + Z 23 ) 2 I(l 12 ,l 2i , hi)I(h3, hi, hi) 



H(hb, he, hc) 1 ^ 2 sin(i(/ a (, + h c + l ac )) 



sin(i(Z bc + l ac - l ab )) 
and 

H(l ab ,l bc ,l ac )= sin( l^+If+k, ) sin ( U+jM^ ) sin ( U=1m±U ) sin( -i ab +hc+iao ) . 
If J12 = j 2 3, l|41|l simplifies to 



s(k) 



r(fc)- 5 /4 7r 



sin(Z 12 )3/ 4 sin(/ 24 ) 1 /2- 

The analogous formula for q = 1 was checked numerically. 

We do not know whether there exist configurations of the type described in 
Figure 2(c) in the spherical case with edge lengths commensurable to and less than 
7r. Neither to do we have a uniform formula, similar to (|36|l . which describes the 
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asymptotics of the 6j symbols in the neighborhood of such a degeneration. The 
case ji2 + J23 + ji3 = r — 2 seems even more mysterious. 

7. COLINEAR VERTICES OR DEGENERATE EDGES 

Suppose that every face of a tetrahedron r is degenerate. We may assume 
without loss of generality that the vertices lie on a geodesic with length at most 7r 
in the order 1, 2, 3, 4 so that j a t + jt c = jac for any l<a<6<c<4. In this case 
the 6j symbol equals 

(_l)3«"([2j 13 + l][2j 24 + IT 1 ' 2 ~ (-l) 2fc ^7rr(fc)- 1 (sina 13 )sm(/ 24 ))- 1 / 2 . 

In the case that one of the lengths is zero one obtains from © the formula 

fci34 ^2 = (- 1 )' £J12 + fej23 + fejl3 ([ 2fc ^ + l][2^23 + l])- 1/2 

_ (_ 1 )fe(ji2+j2 3 +ii3) 7rr ( fc )-i s in(Z 12 )-2 sin^s)^ 

if both lengths are non-zero. Similarly, one obtains the last two formulas in Theo- 
rem 

8. The classically forbidden case 

Naturally one expects 6j symbols are exponentially decreasing in k if r does 
not exist; this was proved that in the case q = 1. The difference equation of Section 
4 together with the results of Geronimo, Bruno and Assche ^21 prove that the 6j 
symbols are exponentially decreasing as a function of a label j a {, after it passes into 
the classically forbidden regime. However, it is not clear to us whether this implies 
exponential decay as a function of k for arbitrary classically forbidden tetrahedra. 

9. The hyperbolic case 

In the case q — exp(7r/r) (no factor of i) the asymptotics are related to 
hyperbolic tetrahedra. If it exists, let r(fc) denote the hyperbolic tetrahedron with 
edge lengths l a b{k) — 2i:(kjab + |)A"(^) an d exterior dihedral angles 6 a b(k). Let 
r denote the hyperbolic tetrahedron with edge lengths l a b = 2n j a ^j [r — 2) and 
exterior dihedral angles 9 a t- One has identities similar to the spherical case, e.g. 

smh(Z cd ) 

where A a — \ det(— cosh(4 c |f )jC7 i a ))| 1 / 2 . We state the results in the hyperbolic case 
as follows. They are identical to the spherical case, except for the replacements sin 
with sinh, cos with — cosh, and Vol(r) with — Vol(r)). The proofs are similar. 



Theorem 9.0.6. Let j t>, l<a<6<4 be non-negative half-integers, r > 2 and 

kj\2 kj 23 kj 13 

kj 3 4 kj U kj 2i , q=cxp(7r/r{k)) 



S(k) 



:!() 
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(a) If t exists and is non- degenerate, then 

2?r cos(</)(fc) + tt/4) 



s(k) 



r(fc) 3 / 2 |det(-cosh(U))| 1 /4 
where 

m = ^ (EM*0U(fc) + 2Vol(r(fc))) . 

n \a<b J 

(b) If t exists, has zero volume but all faces have non-zero area then 

(d) J/ r exists and has exactly one edge length, say l a b vanishing, then 

s{k) ~ (-l) fc ^ + ^ +j! '^7rr(fc)- 1 (sinh(L c )sinh(? M ))- 1 / 2 . 

(e) // r exists, all faces have zero area but all edge lengths are non-vanishing 
then t lies on a geodesic, say with vertices in order a, b, c, d then 

s(k) ~ (-l) 2fc ^^r(fc)- 1 (sinhG ac )sinh(/ M ))- 1 /2. 

(f ) // r exists, all faces are degenerate, and all non-zero edge lengths are equal 
then ( supposing without loss of generality that l ac i =/= Q) 

(g) If t exists, but all edge lengths are vanishing then s(k) = 1. 

10. Questions 

(a) Is there a uniform formula which includes all cases, just as the Schulten- 
Gordon formula (|36|l includes cases (a) and (b)? Such a formula (or least 
uniform estimates) would be needed to apply these results to the asymp- 
totics of the quantum invariants such as Turaev-Viro and Jones, since 
one would presumably have to show that the contributions from the non- 
degenerate case are dominant. 

(b) Is there a geometric description of the tensor category of representations of 
U q (s[2) in terms of the moduli space of hyperbolic bundles [2], for q positive 
real, similar to the description for q a primitive root of unity using unitary 
bundles, which gives a geometric explanation of the hyperbolic formulas? 

(c) Are there results for other values of q, for instance, q negative real, or 
q = exp(7ri/r) but j a b [0, (r — 2)/2]? Numerical experiments show that 
in some of these cases the 6j symbols are rapidly increasing. 

(d) Are there similar formulas for 6j symbols related to the Kashaev-Reshetikhin 
invariants |2l)j? 
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(e) There are analogs of the 6j symbols for other groups. For example, for 
any four-tuple of dominant weights j'12, J23, J34, ju such that the tensor 
products Vj ab ® Vj bc are multiplicity-free, one has two canonical bases for 
the space of invariants Vj 12 ® Vj 23 ® Vj 34 (g) Vj 14 given by the two ways of 
pairing. Are there explicit asymptotic formulas for these symbols? 

Appendix A. Maple code 
The following code (not very well optimized) was used to generate Figure 3. 

with(plots) ; with(linalg) ; 

# quantum integer factorial 

pi :=evalf (Pi) ; nq: =(r ,n) ->sin(pi*n/r)/sin(pi/r) ; 
f aq: =(r ,n) ->evalf (product (nq(r , ii) , ii=l . .n) ,30) ; 

# quantum triangle symbol, Racah sum, quantum sixj symbol 
trq: = (r,a,b,c)->faq(r, (b+c-a) IT) *f aq(r , (a+c-b)/2)* 

faq(r, (a+b-c)/2)/(f aq(r, (a+b+c)/2+l) ) ; 
syq:=(r,al,a2,a3,a4,bl,b2,b3)->sum((-l) ~z*faq(r ,z+l)/ 

(faq(r ,bl-z)*faq(r ,b2-z)*faq(r ,b3-z) *f aq(r ,z-al) * 

faq(r,z-a2)*faq(r,z-a3)*faq(r,z-a4) ) , 

z=max(al,a2,a3,a4) . .min(bl ,b2 ,b3) ) ; 
s j q : = (r , a , b , c , d , e , f ) -> ' if ' (min (a+b-c , a-b+c , -a+b+c , a+e-f , e-a+f , -e+a+f , 

d+c-e , d-c+e , -d+c+e ,b+f -d, b-f +d , -b+f +d , 2*r-4-a-b-c , 2*r-a-e-f , 2*r-d-c-e , 

2*r-b-d-f )>= O.evalf ((trq(r,a,b,c)*trq(r,a,e,f ) *trq(r ,d, c ,e) * 

trq(r,b,d,f))~(l/2)*syq(r, (a+b+c)/2, (a+e+f )/2, (c+d+e)/2, (b+d+f )/2, 

(a+b+d+e)/2, (a+c+d+f)/2, (b+c+e+f )/2) ,30) ,0) ; 

# length corresponding to a dominant weight (integer) 
len:=j->pi*(j+l)/r ; 

# The Gram matrix, its determinant and inverse 

G:=(i, j ,k,l,m,n)->linalg [matrix] (4,4, [1, cos(i) ,cos(k) ,cos(m) , 

cos(i),l, cos(j) ,cos(n) , 

cos (k) , cos ( j ) , 1 , cos(l), 
cos (m) , cos (n) , cos (1) , 1] ) ; 

detG:=(i,j ,k,l,m,n)->evalf (det(G(i, j ,k,l,m,n)) ,10) ; 

Ginv:=(i, j ,k,l,m,n)->evalf (inverse(G(i , j ,k,l,m,n))) ; 

# Does the tetrahedron exist? 

tetexist : = (r , a , b , c , d , e , f ) -> ' if ' (min (a+b-c+1 , a-b+c+1 , -a+b+c+1 , a+e-f +1 , 
e-a+f +1 , -e+a+f +1 , d+c-e+1 , d-c+e+1 , -d+c+e+1 ,b+f -d+1 ,b-f +d+l , -b+f +d+l , 
2*r-4-a-b-c-3 , 2*r-a-e-f -3 , 2*r-d-c-e-3 , 2*r-b-d-f -3 , 
detG(len(a) ,len(b) ,len(c) ,len(d) ,len(e) ,len(f ))) >= 0,1,0); 

# The amplitude of the asymptotic formula 

amppredict : =(r , i , j ,k,l ,m,n)->evalf ( ( (r/pi) "3*1 . 5*pi* (abs (detG( 

len(i) ,len(j),len(k) ,len(l) , len(m) , len(n) ) ) / (36) ) " (1/2) ) ~ (-1/2) , 10) ; 

# The Dihedral angles 

preangle : =(i , j ,k,l,m,n)->evalf (seq(seq(pi - 

arccos(-Ginv(i, j ,k,l,m,n) [a,b]/ (Ginv(i, j ,k,l,m,n) [a, a]* 
Ginv(i,j,k,l,m,n) [b,b] ) " (1/2) ) ,b=a+l . . 4) , a=l . . 4) , 10) ; 
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diangle :=(r , i , j ,k,l,m,n)->evalf (preangle (len(i) ,len( j ) , len(k) , len(l) , 
len(m) , len(n) ) ) ; 

# The "Dehn invariant" aka Regge action 

dehn: =(r , i , j ,k,l,m,n)-> diangle (r ,i,j,k,l,m,n)[l]* (1+1) /2+ 

diangle(r,i, j ,k,l,m,n) [2] *(n+l)/2+ diangle (r , i , j ,k,l,m,n) [3]*(j+l)/2+ 
diangle(r,i, j ,k,l,m,n) [4] *(m+l)/2+ diangle (r , i , j ,k,l,m,n) [5]*(k+l)/2+ 
diangle (r, i, j ,k,l,m,n) [6] *(i+l)/2; 

# The asymptotic formula for the degenerate case 
G3: = (i, j ,k)->linalg [matrix] (3,3, [l,cos(i) ,cos(j) , 

cos (i) , 1 , cos (k) , 
cos ( j ) , cos (k) , 1] ) ; 
areaq:=(r,i, j ,k)->evalf (det(G3(len(i) ,len(j) ,len(k) ) ) ) " (1/2) ; 
tang:=(r,i, j ,k,l ,m,n)->evalf (r~ (-4/3) *2~ (2/3) *3~ (-1/3) *pi~ (4/3) * 
GAMMA (2/3) " (-1) * (areaq(r , i , j ,k) *areaq(r , i ,m,n) *areaq(r ,1 , j ,n)* 
areaq(r,l,m,k))~(-l/6)) ; 

# NUMERICAL EXPERIMENT 

r:=200; i:=40; j:=48; k:=50; 1:=52; m:=54; stepsize : = . 1 ; 
u: =2*max(i , j ,k,l ,m) ; 
#calculate the 6j symbols 

symbols:= [ seq( [2*y , sjq(r , i , j ,k, 1 ,m, 2*y)] ,y=0 . .u/2) ]; 
#does the terahedron exist 

exist : = [seq( [y*stepsize ,tet exist (r , i , j ,k, 1 ,m,y*stepsize)] , 

y=l. . (r/2)/stepsize)] ; 
#the predictions for the tangent case, plus or minus 
tapredict:= [ seq( [y , ' if ' (tetexist (r , i , j ,k, 1 ,m,y)=l , 

tang(r,i, j ,k,l,m,y) , undefined)] ,y=0. .u/2) ] ; 
mtapredict:= [ seq( [y, -' if ' (tetexist (r, i , j ,k, 1 ,m,y)=l , 

tang(r,i, j ,k,l,m,y) , undefined)] ,y=(u/2) . .u) ] ; 
#the dihedral angles 

angpredict := [ seq( [y*stepsize .diangle (r, i , j ,k,l ,m, 

y*stepsize) [2]] ,y=0. .u/stepsize) ] ; 
for y from to u/stepsize do 

phqpredict [y+1] := ' if ' (tetexist (r, i ,j ,k,l ,m,y*stepsize) = l , 
phqpredict [y] + angpredict [y+1] [2] *stepsize/2 , 
'if (tetexist (r,i, j ,k,l,m, (y+1) *stepsize) = l , 
dehn(r , i , j ,k,l ,m, (y+1) *stepsize) ,0) ) od; 
#the predicted amplitudes , plus or minus 

predict a: = [ seq( [y*stepsize , 'if (tetexist (r , i , j ,k, 1 ,m,y*stepsize)=l , 

amppredict(r,i, j ,k,l,m,y*stepsize) .undefined)] ,y=l. .u/stepsize) ] ; 

mpredicta: = [ seq( [y*stepsize , ' if ' (tetexist (r , i , j ,k, 1 ,m,y*stepsize)=l , 

-amppr edict (r , i , j ,k, 1 ,m,y*stepsize) .undefined)] ,y=l . .u/ stepsize) ] ; 

#the prediction for the non-degenerate case 

predict:= [seq( [stepsize*y ,predicta[y] [2] *cos (pi* (1/4) + 
phqpredict [y] )] ,y=l . .u/stepsize)] ; 

# Show all the plots 

plot ( [symbols , predict ,mtapredict , tapredict .predict a, mpredicta] , 
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style= [point , line , line , line , line , line] , 

# color= [red , blue , green , green , yellow , yellow] , 

# linestyle=[SOLID ) SOLID ) DOT,DOT > SOLID,SOLID] , 
symbol=circle , symbolsize=5) ; 
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